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ABSTRACT: The Cassini Division in Saturn's rings contains a series of eight named 
gaps, three of which contain dense ringlets. Observations of stellar occultations by the 
Visual and Infrared Mapping Spectrometer onboard the Cassini spacecraft have yielded 
~40 accurate and precise measurements of the radial position of the edges of all of these 
gaps and ringlets. These data reveal suggestive patterns in the shapes of many of the gap 
edges: the outer edges of the 5 gaps without ringlets are circular to within 1 km, while 
the inner edges of 6 of the gaps are eccentric, with apsidal precession rates consistent with 
those expected for eccentric orbits near each edge. Intriguingly, the pattern speeds of these 
eccentric inner gap edges, together with that of the eccentric Huygens ringlet, form a series 
with a characteristic spacing of 0.06°/day. 

The two gaps with non-eccentric inner edges lie near first-order Inner Lindblad Reso- 
nances (ILRs) with moons. One such edge is close to the 5:4 ILR with Prometheus, and 
the radial excursions of this edge do appear to have an m = 5 component aligned with that 
moon. The other resonantly confined edge is the outer edge of the B ring, which lies near the 
2:1 Mimas ILR. Detailed investigation of the B-ring-edge data confirm the presence of an 
m = 2 perturbation on the B-ring edge, but also show that during the course of the Cassini 
Mission, this pattern has drifted backward relative to Mimas. Comparisons with earlier 
occultation measurements going back to Voyager suggest the possibility that the m = 2 
pattern is actually librating relative to Mimas with a libration frequency L ~ 0.06°/day 
(or possibly 0.12°/day). In addition to the m = 2 pattern, the B-ring edge also has an 
m = 1 component that rotates around the planet at a rate close to the expected apsi- 
dal precession rate {vjb ~ 5.06°/day). Thus the pattern speeds of the eccentric edges in 
the Cassini Division can be generated from various combinations of the pattern speeds of 
structures observed on the edge of the B ring: £l p = tdb — jL for j = 1,2,3, 7. We 
therefore suggest that most of the gaps in the Cassini Division are produced by resonances 
involving perturbations from the massive edge of the B ring. We find that a combination 
of gravitational perturbations generated by the radial excursions in the B-ring edge and 
the gravitational perturbations from the Mimas 2:1 Inner Lindblad Resonance yields terms 
in the equations of motion that should act to constrain the pericenter location of particle 
orbits in the vicinity of each of the eccentric inner gap edges in the Cassini Division. This 
alignment of pericenters could be responsible for forming the Cassini-Division gaps and thus 
explain why these gaps are located where they are. 

Subject Keywords: occultations, planets: rings 
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Introduction 



The Cassini Division is a roughly 4500-km wide region in Saturn's main rings situated 
between the A and B rings. Far from being a completely empty gap between these two 
rings, this zone is actually a complex region containing an array of gaps and ringlets. The 
physical processes responsible for creating and maintaining most of the observed features 
in this region remain obscure. 

In particular, there is still no definitive explanation for most of the numerous gaps 
present throughout the inner part of Cassini Division. As shown in Figure [H there are eight 
gaps in the inner part of the Cassini Divisi on, all of which are now named after various 
researchers who worked on Saturn's rings (jColwell et al.l l2009bl ). The innermost gap is 
called the Huygens Gap and it marks the inner boundary of the Cassini Division. The 
inner edge of the Huygens Gap, which is also the outer edge of the massive B ring, has long 
been known to be associated with a 2:1 mean-motion resonance with Saturn's moon Mimas 



(jPorco et al.lll984l ). This gap contains two optically thick ringlets: the inner one, called the 
"Huygens Ringlet" , is known to be eccentric, while the outer, or "Strange" ringlet seems t o 
have a significant inclination (jTurtle et al.lll99ll : ISpitale fc Porcdl2006l : ISpitale et al.ll2008l ). 



While the Mimas 2:1 resonance likely plays an important role in creating the Huygens 
Gap, the origins of the other seven gaps, as well as the dense ringlets in the Huygens, 
Herschel and Laplace gaps, are much l ess clear. Some hav e argued that each gap in the 
Cassini Division contains a tiny moon (jLissauer et al.lll98ll ). just as the Encke and Keeler 
Gaps in the A ring are maintained by the small moons Pan and Daphnis. While some 
wavelike featur es in the Cassini Divisi on have been interpreted as evidence for the existence 
of such moons (jMarouf Tvlerlll986l ). direct detections of the moons themselves have not 
yet been reported. It is also not clear why such moons would be concentrated in this 
particular part of the ring system, well inside Saturn's Roche limit for icy bodies. 

Using the extensive occultation data obtained by the Visual and Infrared Mapping 
Spectrometer (VIMS) onboard the Cassini spacecraft, we have conducted an investigation of 
the gaps in the Cassini Division. The high resolution and accurate geometrical information 
possible with these data has enabled us to determine the shapes of most of the edges of 
these gaps and ringlets. Based on this information, we will propose alternative explanations 
for these gaps that do not require the existence of a small moon within each and every gap. 

Section 2 of this paper described the VIMS observations used in this analysis, and how 
they were processed to obtain position estimates for the various edges. Section 3 presents 
the results of our analysis and the derived constraints on the shapes of various edges in the 
Cassini Division, including the B-ring edge. The B-ring edge turns out to be quite complex, 
so Section 4 briefly compares some of this edge's observed features to the predictions from 
simple dynamical models. Section 5 presents a dynamical model that could explain the 
location of the inner edges of the various gaps in the Cassini Divsion as the result of a series 
of resonances involving perturbations from the massive B ring edge. This model is only 
a first step towards understanding the architecture of the Cassini Division, so Section 6 
discusses potential future work that could confirm and extend this model. The final section 
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Fig. 1. — Overview of the basic architecture of the inner part of the Cassini Division. This 

optical depth profile is derived from the Rev 8 ingress occultation of o Ceti by the rings 

observed by the VIMS instrument on board Cassini. The various gaps discussed in this 
study are labeled. 

summarizes our conclusions. 



2. Observations and Data Reduction 

The Visual and Infrared Mapping Spectrometer (VIMS) is most often used to produce 
spatially-resolved spectra of planetary targets between 0.3 and 5.1 microns. However, VIMS 
is a fle xible instrument t hat can operate in a variety of modes, including an occultation 



mode (IBrown et al.1 12004). In this mode, the imaging capabilities are disabled, the short- 
wavelength VIS channel of the instrument is turned off and the IR channel obtains a series 
of 0.8-5.1 micron spectra from a single pixel targeted at a star. Typical sampling invtervals 
were 20-80 ms. The data used in this analysis are uncalibrated, but a mean instrumental 
thermal background has been subtracted from all the spectra for each occultation. A precise 
time stamp is appended to every spectrum to facilitate reconstruction of the occultation 
geometry. 



2.1. Observations 



Up through the end of 2008, VIMS has observed over 50 occultations, which have 
yielded a total of 48 potentially useful cuts through the Cassini Division. Table [U lists the 
occultation cuts used in this analysis, along with the occultation times and inertial longi- 
tudes of the observations, the maximum Data Number (DN) detected, the radial sampling 
scale and the shift required to bring circular features into alignment (see below). 
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Table 1: Occu pation cuts used in this analysis 



Rev 


Star 




UTC 


Inert ial 


max DX 


Radial 


Radial 


Notes 








at Jeffreys Gap 


Longitude 


at 3 /iin 


Sampling (km) 


Offset (km) 




08 


oCet 


1 


2005-144T05:58 


1° 


995 


0.30 


+0.06 


A 


08 


oCet 




2005-144T07:07 


-28° 


999 


0.30 


-0.68 


A 


13 


qSco 


i 


2005-232Tll:46 


-70° 


860 


0.37 


+0.77 


B.I 


13 


— 




2005-232T13:44 


-3° 


770 


0.38 


-0.98 


I 


26 






2005-204T16:47 


-3° 


967 


0.21 


-0.32 




28 


aTau 


i 


2006-252T10:54 


34° 


144 


0.52 


+0.15 


F 


29 


<SVir 


1 


2006-268T22:37 


204° 


125 


1.44 


-1.35 


Xr,A,C 


29 


<SVir 




2006-268T22:44 


97° 


128 


1.43 


+2.74 


Xr,A,C 


29 




. 


2006-269T07:41 


-166° 


712 


0.10 


+0.24 




30 


RLeo 


* 


2006-285T02:ll 


-39° 


75 


0.44 


+0.59 


Xd,C,I 


30 


RLeo 




2006-285T02:39 


-85° 


74 


0.44 


-0.08 


Xd,C,I 


31 


CWLeo 


i 


2006-301T01:26 


-23° 


167 


3.07 


+1.21 


Xr.B 


31 


CWLeo 


e 


2006-301T02:10 


-103° 


188 


3.07 


-0.91 


Xr.B 


34 


ojAur 




2006-336T12:54 


27° 


414 


0.92 


+0.21 


F 


36 


RHya 


j 


2007-001T17:07 


-165° 


333 


0.35 


+0.43 


B 


41 




j 


2007-082T17:30 


8° 


212 


0.31 


-0.45 


B,I 


41 


RHya 


j 


2007-088T07:24 


-156° 


112 


0.09 


+0.18 


B,D,I 


63 


RLeo 




2008-094T13:ll 


96° 


368 


0.23 


+0.93 




63 


RLeo 




2008-094T13:58 


126° 


371 


0.23 


+0.01 


G 


65 


RCas 




2008-112T00:43 


30° 


77 


0.30 


+0.27 


Xd,C,I 


70 


CWLeo 




2008-155T14:22 


82° 


395 


0.76 


-2.11 




70 


CWLeo 




2008-155T16:15 


137° 


363 


0.75 


-4.15 




71 






2008-160T09:04 


-173° 


422 


0.28 


+0.46 




72 


7 Cru 




2008-167T12:21 


-174° 


639 


0.28 


+0.83 




73 


7Cru 




2008-174T15:30 


-174° 


635 


0.28 


+0.46 




75 


RLeo 




2008-191T04:43 


80° 


279 


0.36 


+1.07 




75 


RLeo 




2008-191T06:27 


138° 


279 


0.36 


-0.63 




77 


RLeo 




2008-205T06:58 


84° 


295 


0.31 


+0.10 




77 


RLeo 




2008-205T08:27 


133° 


296 


0.31 


+0.18 




78 






2008-209T20:07 


-175° 


304 


0.14 


+0.52 


B 


78 


/3Gru 




2008-210T09:35 


-96° 


254 


0.34 


-0.26 




80 


RSCnc 


* 


2008-226T02:03 


57° 


310 


0.56 


-0.61 


H 


80 


RSCnc 




2008-226T07:29 


155° 


305 


0.56 


+0.00 




81 






2008-231T06:57 


-177° 


579 


0.27 


+0.12 




82 


-yCm 


1 


2008-238T15:32 


-177° 


721 


0.27 


-0.04 




85 


RSCnc 




2008-262T22:29 


59° 


306 


0.54 


-0.30 




85 


RSCnc 


e 


2008-263T03:48 


153° 


312 


0.54 


+0.31 




86 


7Cru 


i 


2008-268T03:12 


-178° 


1020 


0.41 


-0.33 




87 


RSCnc 


i 


2008-277T16:16 


59° 


326 


0.53 


-0.12 




87 


RSCnc 


e 


2008-277T2L30 


152° 


325 


0.53 


+0.01 




89 


7Cru 




2008-290T04:23 


-178° 


704 


0.27 


-0.13 


B 


92 


RSCnc 




2008-315T0L56 


85° 


225 


0.21 


+0.67 


B 


93 


7Cru 




2008-320T16:31 


-157° 


484 


0.18 


+0.09 


B.E 


94 


7Cru 




2008-328T01:23 


-168° 


267 


0.13 


+0.30 


B.E 


94 


fMus 




2008-328T08:20 


-102° 


188 


0.19 


+0.53 


B 


94 


eMus 


e 


2008-328T13:02 


-44° 


188 


0.19 


+0.02 


B 


96 


7Cru 


i 


2008-343Tll:45 


-172° 


236 


0.14 


+0.18 


B 


97 


7Cru 


i 


2008-351Tll:03 


-172° 


794 


0.42 


+0.44 





Note times, longitudes and maximum DN are evaluated inside the Cassini Division. 
Xr: Excluded from Cassini Division analysis due to coarse radial sampling 
Xd: Excluded from Cassini Division analysis due to low signal levels 

A: Low inclination occs, Saturn pole position adjusted to achieve approximate match with circular 
features 

B: Data smoothed by 3 prior to edge-finding 

C: Data resampled and smoothed by 5 prior to edge-finding 

D: Data smoothed by 10 prior to finding Bcssel Gap Inner Edge 

E: Data smoothed by 12 prior to finding Laplace Gap Inner Edge 

F: Jeffreys Gap Inner edge falls in data gap 

G: Laplace Ringlet Inner edge falls in data gap 

H: Kuiper Gap Inner edge falls in data gap 

I: Spectral channels covering the range from 2.73-3.11 microns used. 
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Table 2: Assumed star positions in the International Celestial Reference Frame at the Hip- 
parcos epoch. Data obtained from the Simbad astronomical database and the Hippacos 
online catalog. Note that all the stars except CWLeo and ?7Car were present in the Hip- 
parcos catalog. 



Star 


Right Ascension 


Declination 


Proper Motion 


Proper Motion 


Parallax" 




hr:min:sec 


hr:min:sec 


R.A. (mas/yr) 


Dec. (mas/yr) 


mas 


oCet 


02:19:20.7866 


-02:58:37.418 


+10.33 


-239.48 


7.79 


aSco 


16:29:24.4675 


-26:25:55.006 


-10.16 


-23.21 


5.40 


aOri 


05:55:10.2892 


+7:24:25.331 


+27.33 


+10.86 


7.63 


aTau 


04:35:55.2005 


+16:30:35.142 


+62.78 


-189.36 


50.09 


<5Vir 


12:55:36.4833 


+03:23:51.355 


-471.44 


-52.81 


16.11 


RLeo 


09:47:33.4907 


+11:25:44.020 


-0.57 


-42.70 


9.87 


CWLeo 


09:47:57.382 


+13:16:43.66 








aAur 


05:16:41.2956 


+45:59:56.505 


-75.52 


-427.13 


77.29 


RHya 


13:29:42.8189 


-23:16:52.888 


-60.73 


+11.01 


1.62 


RCas 


23:58:24.7936 


+51:23:19.545 


+84.39 


+18.07 


9.37 


7Cru 


12:31:09.9293 


-57:06:45.249 


+27.94 


-264.33 


37.09 


r/Car 


10:45:03.591 


-59:41:04.26 


-7.6 


+1.0 




/3Gru 


22:42:39.9349 


-46:53:04.437 


+ 135.68 


-4.51 


19.17 


RSCnc 


9:10:38.8054 


+30:57:47.589 


-9.41 


-33.05 


8.21 


eMus 


12:17:34.6363 


-67:57:38.418 


-231.26 


-26.37 


10.80 



a Parallax measured at Earth 
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2.2. Geometrical Navigation 

Table[2]gives the assumed positions of the stars used in this analysis in the International 
Celestial Reference Frame. For each occultation, the position of the star is adjusted to 
account for both the proper motion of the star and the parallax at Saturn. The available 
SPICE kernels provided by the Cassini navigation team were then used to predict the 
apparent position (radius and inertial longitude) of the star in Saturn's ring plane as a 
function of time in a planetocentric reference frame, taking into account stellar aberration. 
In nearly all cases, this estimate of the occultation geometry was confirmed to be accurate 
to within a few kilometers using the known radii of nearly circular gap edges in the Cassini 
Division and the outer A Ring (Features 1, 3, 4, 13, 16 and 20 of French et al. 1993). 
The exceptions were the low- inclination stars o Ceti and 5 Virginis, for which features 
could be tens of kilometers away from their nominal positions. In these cases, the fiducial 
position of Saturn's pole was adjusted slightly (by at most 0.015°) to bring these cuts into 
alignment with the other occultations. The residual scatter in the estimated locations of 
circular features is consistent with uncertainties in the star positions, observation timing 
and spacecraft trajectory (see below). 

2.3. Lightcurve Generation 

While each occultation produces several time series of brightness measurements at 
multiple wavelengths, only the average brightness around a wavelength of three microns 
is used for this analysis. This is because in general, the signal measured by the VIMS 
instrument was a combination of transmitted starlight and reflected sunlight. The reflected 
signal from the rings is strongly attenuated at wavelengths close to the strong 2.9 micron 
water-ice absorption band, so using the spectral channels in this range minimizes the ring 
background. 

To save on data volume, in many of these occultations the normal spectral resolution 
of the instrument was reduced by a factor of eight. A total of 32 spectral channels were 
returned with an average resolution AA ~ 0.13/xm. For the few occultations taken at 
full spectral resolution, the data from the appropriate spectral channels were co-added 
in software after the fact to make their spectral resolution and signal-to-noise consistent 
with the other occulations. For most of the occultations, the summed spectral channel 
corresponding to the wavelength range between 2.87 and 2.98 microns was used in this 
analysis. However, a small number of occultations (Rev 13 a Scorpii, Rev 30 R Leonis, 
Rev 41 a Auriga, Rev 41 R Hydra and Rev 65 R Cassiopaea), had either sufficiently low 
signal-to-noise ratios or sufficiently high cosmic-ray-fluxes that a single profile was too noisy 
to obtain reliable edge detections. In these cases, three brightness profiles were averaged 
together, covering the wavelength range of 2.73-3.11 microns. This simple averaging reduced 
the instrumental noise and cosmic ray background to acceptable levels. 

While using data near 2.9 microns minimizes the background from the illuminated rings, 
it does not completely eliminate this signal from all occultations. This residual ring back- 
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groun d can be eliminated by comparing data at different wavelengths ([Nicholson fc Hedman 



20091 ) , but such refinements are not necessary here because the goal of this analysis is simply 



to determine the locations of sharp edges. The time variations in the ring background are 
limited by the spatial resolution of the instrument, while the time variations in the stellar 
signal are limited by the temporal resolution, so the ring background varies slowly with 
time compared to the variations in the stellar signal and has negligible effect on the inferred 
locations of sharp edges. 



2.4. Edge Detection 

After deriving the lightcurve and geometrically navigating each occultation, we deter- 
mined the radius, inertial longitude and time when the star passed behind 22 edges in the 
Cassini Division: the inner and outer edges of the Huygens, Herschel, Russell, Jeffreys, 
Kuiper, Laplace, Bessel and Barnard Gaps, as well as the inner and outer edges of the 
dense ringlets in the Huygens, Herschel and Laplace Gaps. Note that we do not consider 
the "Strange" ringlet in the Huygens Gap in this analysis, nor the low-optical depth ringlets 
in the Huygens, Jeffreys or Laplace Gaps. 

The edges in the Cassini Division are quite variable in their morphology, and even 
a single edge may have different shapes in diff erent occultations . This could make the 



identification of edges based on half-light levels ([French et al.lll993l ) problematic, so instead 
we identified the edge as the point in the lightcurve with the steepest slope. Even so, 
some effort was required to prevent the edge-detection algorithm from mistaking cosmic ray 
spikes or other sharp ring features for the desired edge. The algorithm therefore begins by 
normalizing the total signal (star + reflected ringshine, if any) to unity in the gaps. This is 
accomplished by first determining the average DN in two clear zones (117,700-117,750 km 
and 119,970-120,000 km) within the two widest gaps in the Cassini Division: the Huygens 
Gap and the Laplace Gap. A linear trend based on these two numbers establishes the 
unocculted stellar brightness throughout this region and is used to normalize the data so 
that the signal in each gap is approximately 1.0. 

For each edge, a fiducial zone is selected for analysis which is sufficiently wide to 
accommodate all observed edge positions and residual pointing errors. (For the B ring edge 
this region is 150 km wide. For the Huygens ringlet edges this zone is 40 km wide. For the 
Herschel Ringlet and Herschel Gap inner edge it is 20 km wide. For the inner edge of the 
Laplace gap it is 15 km wide. For the Herschel gap outer edge, the Russell and Jeffreys 
gap inner edges, the Laplace ringlet and the Laplace gap outer edge, and the Bessel and 
Barnard gaps it is 10 km. For the Russell and Jeffreys gaps' outer edges and the Kuiper 
gap edges it is 5 km.) The algorithm first makes a preliminary estimate of the location of 
the relevant gap edge based on where the signal in this region first deviates significantly 
from unity. The final estimate of the edge location is then the point where the brightness 
profile has the steepest slope in a region within ±4 km of the preliminary edge estimate 
(except for the low-resolution 5 Virginis occultation, where a region ±20 km was used). 



- 9 - 



For certain occultations with fast time-sampling or low signal-to-noise, the raw profile 
is too noisy for the above algorithm to find all the edges reliably. In these cases, a boxcar 
smoothing was applied to the data prior to estimating the edge position. The occultations 
where this was done and the specific smoothing lengths used are given in Table [TJ In all 
cases, the uncertainty in the edge position is determined by the radial sampling interval 
because this is typically much larger than the projected stellar diameter or the fresnel zone. 

2.5. Occultation Quality 

While the above procedure was applied to all the occultation cuts in Table [fl there 
were a few occultations which we elected not to use with the Cassini Division edges. The 
Rev 30 R Leonis and Rev 65 R Cassiopaea occultations had signal-to-noise ratios too low 
to detect all the edges in the Cassini Division reliably, while the resolution of the Rev 29 5 
Virginis and the Rev 31 CW Leonis occultations were so much lower than the rest of the 
observations that they would not contribute much to our understanding of most of these 
edges. These data are therefore not considered "Quality" occultations and are not used 
in the analysis involving the edges within the Cassini Division itself (they are designated 
with Xr or Xd in Table [1]). However, many of these occultations occur at early times in the 
Cassini mission, when VIMS occultations are comparatively rare, so they are useful when 
considering the time-evolution of the B ring edge (see below) . 

3. Data Analysis 

The above procedures yield roughly 40 "quality" estimates of the position of each of 
the 22 selected gap edges for a range of different times and inertial longitudes. The simplest 
statistics that can be computed from these data are the (unweighted) means and standard 
deviations of all the edge estimates, which are presented in Table El Of particular interest 
are the standard deviations, which range from around 1 km up to 50 km for the edge of the 
B ring. 

Based on these variances, we can begin to identify several different classes of edges 
in the Cassini Division. First of all, there are several features whose scatter is around 1 
km, which likely represent truly circular ring edges. These include the outer edges of the 
Russell, Jeffreys, Kuiper, Barnard and possibly the Bessel and Herschel Gaps. (The inner 
edge of the Kuiper gap has a similarly low variance, but further investigation shows that it 
belongs in the eccentric class of ring features, discussed below). 

The remaining edges all have noticeably larger variances and therefore appear to be 
non-circular in some form. Three of these edges (the outer edge of the B ring, and the inner 
edges of the Laplace Ringlet and the Barnard Gap) are close to mean-motion resonances 
with known satellites of Saturn (Mimas 2:1, Pandora 9:7 and Prometheus 5:4, respectively) 
which could provide natural explanations for their shapes. However, this still leaves most of 
the edges unexplained. We have found that 10 of these features (both edges of the Huygens 



-10- 



Table 3: Elementary Properties of Cassini Division Gap Edges 



Feature 


Mean Radius 


St.Dev. 


Edge Type 




(km) 


(km) 




B-ring Outer Edge 


117564.4 


51.4 


Resonant (2:1 Mimas) 


Huygens Ringlet Inner Edge 


117804.4 


20.3 


Eccentric 


Huygens Ringlet Outer Edge 


117822.8 


20.2 


Eccentric 


Huygens Gap Outer Edge 


117930.6 


2.8 


Unknown 


Herschel Gap Inner Edge 


118188.2 


5.9 


Eccentric 


Herschel Ringlet Inner Edge 


118233.9 


2.3 


Unknown 


Herschel Ringlet Outer Edge 


118263.2 


2.7 


Unknown 


Herschel Gap Outer Edge 


118283.4 


1.1 


Unknown (Circular?) 


Russell Gap Inner Edge 


118589.7 


5.2 


Eccentric 


Russell Gap Outer Edge 


118628.2 


0.9 


Circular 


Jeffreys Gap Inner Edge 


118929.6 


2.7 


Eccentric 


Jeffreys Gap Outer Edge 


118966.5 


0.8 


Circular 


Kuiper Gap Inner Edge 


119401.7 


1.0 


Eccentric 


Kuiper Gap Outer Edge 


119406.1 


0.7 


Circular 


Laplace Gap Inner Edge 


119845.2 


2.6 


Eccentric 


Laplace Ringlet Inner Edge 


120036.4 


1.8 


Unknown (Resonant? 9:7 Pandora) 


Laplace Ringlet Outer Edge 


120077.7 


2.1 


Eccentric 


Laplace Gap Outer Edge 


120085.7 


1.3 


Eccentric 


Bessel Gap Inner Edge 


120231.3 


1.6 


Eccentric 


Bessel Gap Outer Edge 


120243.6 


1.0 


Circular? 


Barnard Gap Inner Edge 


120303.7 


2.4 


Resonant? (5:4 Prometheus) 
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ringlet, the inner edges of the Herschel, Russell, Jeffreys, Kuiper, Laplace and Bessel Gaps, 
and the outer edges of the Laplace gap and ringlet) can be well fit with ellipses that precess 
around the planet at the expected rate given Saturn's oblateness. Four non-circular edges 
(both sides of the Herschel ringlet, as well as the outer edges of the Huygens and possibly 
the Herschel Gaps) cannot be fit by simple elliptical models and therefore require further 
investigation. 

We will consider each of these different classes in detail below. First, we will investigate 
the apparently circular edges and show how they can be used to refine the estimates of 
the other edge locations. Next, we will discuss the eccentric features and compare the 
amplitudes and precession rates of these features. We will then consider the inner edge of 
the Barnard Gap and explore whether the structure of this edge can be explained by the 
perturbations from Prometheus. Finally, we will look at the outer edge of the B ring, which 
has the largest radial excursions and the most complex structure. The inner edge of the 
Laplace ringlet, both edges of the Herschel ringlet, and the outer edges of the Herschel and 
Huygens Gaps still require further study and will not be discussed in detail in this paper. 



3.1. Circular Edges 

Figure [2] shows the scatter in the inferred positions of the outer edges of the Herschel, 
Russell, Jeffreys, Kuiper, Bessel and Barnard Gaps. All of these data sets have rather 
low dispersions. Furthermore, there is a strong correlation among the radial excursions of 
different edges from the same occultation. (For example, the low points in all six panels 
near 80° and 140° longitude come from the Rev 70 CW Leo occultations. Since this object 
was not observed by Hipparcos, the occultation geometry here is more uncertain than the 
others.) This suggests that most of the scatter in the positions of these edges is due to 
small errors in the geometric reconstructions of different occultations. 

Assuming that most of these edges are circular, we use these features to compute a 
radial offset for each cut and refine the location estimates for the remaining edges. For 
each occultation, we compute the mean deviation of the measured positions of the outer 
edges of the Russell, Jeffreys, Kuiper, Bessel and Barnard gaps from the mean values listed 
in Table [3] (118628.2, 118966.5, 119406.1, 120243.6 and 120315.9 km, respectively). These 
radial shifts, tabulated in Table [H are then applied to all edges for that occultation. Note 
that the nominal positions of the outer edges of the Russell, Jeffreys and Bessel gaps are 
1-2 km interior to the values reported in French et al. (1993). Such a small shift should be 
of no major consequence for this analysis, but will be the subject of a future investigation 
of Saturn's pole orientation incorporating data from all circular features in the rings. 
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Fig. 2. — The radial excursions (deviations from the unweighted mean) of six apparently 
circular edges in the Cassini Division, plotted as functions of inertial longitude. The error 
bars indicate the radial sampling scale of each occultation. 
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3.2. Eccentric Edges 

The observed radial position r' of an eccentric edge depends on the inertial longitude 
A and time relative to some epoch St of the observation as follows: 

r' = r' - A*cos(\-tt p 5t- A ), (1) 

where r' Q is the mean edge location, A is the amplitude of the radial excursions, A Q is a 
phase offset and £l p is a pattern speed. Since an eccentric edge closely follows the path of a 
freely precessing particle on an eccentric orbit, we expect Q p = w, where w is the apsidal 
precession rate, which in the Cassini Division ranges between 4.5°/day and 5.1°/day. 

Preliminary investigations showed that ten of the edges in the Cassini Division could be 
well fit by such an eccentric model with f2 p close to the expected value of w. These features 
are both edges of the Huygens ringlet, the inner edges of the Herschel, Russell, Jeffreys, 
Kuiper, Laplace and Bessel Gaps, as well as the outer edges of the Lapla ce Gap and ringlet . 



The Huygens ringlet w as previously known to be an eccentric structure ([Turtle et al 



1991 



Spitale &: Pored 120061 ) , and there was som e evidence that the inner edge of the Herschel 



Gap might be as well (jFlvnn &; Cuzzilll989l ). However, the fact that a majority of the edges 



in the Cassini Division are simple ellipses is a rather unexpected finding. 

Having identified these features, we sought to determine the parameters in Equation [TJ 
Given the extensive occultation data available, and mindful th at the precise precession r ates 



at a given location might be affected by nearby ring material (INicholson Porcdll988l ). we 
chose not to assume a predicted pattern speed for each edge but instead included Q p as a 
free parameter in each fit. 

For each possible value of the pattern speed fL, we computed the following quantities: 

1 " 

OCR = ~ J2( r 'i ~ f * COS ( A i - Q P St i) ( 2 ) 
1 " 

ui = ~ - f') * sin(Ai - QpSti) (3) 

8=1 

where f' is the mean radius of the edge and the sums are over the ~40 measurements of r'. 
If the radial position of the edge is described by Equation[IJ then J a 2 R + aj = pA/2, where 
p is the correlation coefficient between r\ — r' and cos(Aj — Qpdti — A D ). This function will 
be at a maximum when there is perfect correlation between these parameters (i.e. when 
p = 1). In the limit of perfect sampling of all possible true anomalies (Aj — £l p 5ti), this 
should only occur when the assumed pattern speed Q p equals the true pattern speed. Thus 



the value of Q p that yields the highest value of J a 2 R + aj provides the best estimate of the 
true pattern speed. In this case, the parameters or and a\ will approach — (A/2) cos A 
and —(A/2) sin A (again in the limit of perfect sampling of all possible true anomalies), so 

we can estimate the amplitude and phase parameters using the equations A/2 = a\ + aj 
and tanA = aj/aR. 
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Fig. 3. — Radial excursions (after offset subtraction) of the inner and outer edges of the 
Huygens Ringlet. Error bars indicate the radial sampling scales of the different occultations. 
Note time is measured relative to an epoch of 2005-195T02:13:13.557 (which corresponds 
to a Cassini spacecraft clock time of 150000000), and the mean radii are for the eccentric 
model rather than the data. 
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Fig. 4. — Radial excursions (after offset subtraction) of the inner edges of the Herschel, 
Russell, Jeffreys and Kuiper gaps. Error bars indicate the radial sampling scales of the 
different occultations. Note time is measured relative to an epoch of 2005-195T02: 12: 13.557 
(which corresponds to a Cassini spacecraft clock time of 150000000), and the mean radii 
are for the eccentric model rather than the data. 
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Fig. 5. — Radial excursions (after offset subtraction) of the inner edges of the Laplace and 
Bessel Gaps, as well as the outer edges of the of the Laplace Gap and Ringlet. Error bars 
indicate the radial sampling scales of the different occultations. Note time is measured 
relative to an epoch of 2005-195T02:12:13.557 (which corresponds to a Cassini spacecraft 
clock time of 150000000), and the mean radii are for the eccentric model rather than the 
data. 
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Table 4: Parameters for the Eccentric Edges 



Feature 


Mean Radius" 


Amplitude 


Phase 6 


Pattern Speed 
(Observed) 


Pattern Speed 
(Predicted) c 


Huygens Ringlet Inner Edge 


117805.8 km 


26.9 km 


-80° 


5.03°/day 


5. 022° /day 


Huygens Ringlet Outer Edge 


117823.5 km 


26.9 km 


-74° 


5.03° /day 


5.019° /day 


Herschel Gap Inner Edge 


118188.3 km 


7.9 km 


-171° 


4.97°/day 


4.964° /day 


Russell Gap Inner Edge 


118589.8 km 


7.3 km 


118° 


4.93°/day 


4.904° /day 


Jeffreys Gap Inner Edge 


118929.6 km 


3.4 km 


-98° 


4.87°/day 


4.854° /day 


Kuiper Gap Inner Edge 


119401.6 km 


1.0 km 


13° 


4.79°/day 


4.786° /day 


Laplace Gap Inner Edge 


119844.7 km 


3.6 km 


7° 


4.73°/day 


4.723° /day 


Laplace Ringlet Outer Edge 


120077.9 km 


2.8 km 


123° 


4.72°/day 


4.691°/day 


Laplace Gap Outer Edge 


120085.6 km 


1.5 km 


9° 


4.73°/day 


4.689° /day 


Bessel Gap Inner Edge 


120231.2 km 


2.1 km 


14° 


4.68°/day 


4.669° /day 



a Note these mean radii differ slightly from those in Table [3] because these are the mean 
radii of the eccentric model rather than the mean of the data. 

b Longitude of edge's pericenter for an eccentric model using the given pattern speed at 
an epoch of 2005-195T02:12:13.557 (which corresponds to a Cassini spacecraft time of 
150000000). 

c Predicted precession rate of eccentric particle orbit with semi-major axis equal to the 
observed mean edge r adius, based on current estimates of Saturn's gravity field parameters 
(jjacobson et al.ll2006l ) 



Figures ® [5] show the measured radial excursions for the eccentric edges as a function of 
A — Q p 5t using the best-fit pattern speeds. In all cases, the data are well organized into a sine 
wave. Table 0] lists the best fit parameters for these edges. Note that the observed pattern 
speeds of these features are not far from the predicted pattern speeds derived by assuming 
the edge follows a freely precessing elliptical orbit. The uncertainty in the observed pattern 
speeds is set by the time baseline of the observations, which is approximately 1300 days. A 
difference of 0.01°/day in the pattern speed would therefore shift the phase of the first and 
last occultations by roughly 15°, which is probably on the edge of detectability. 

What is particularly interesting about the fitted pattern speeds of the eccentric edges 
is that they almost form a regular sequence. To see this, first note that both edges of 
the Huygens Ringlet have the same pattern speed (5.03° /day), as expected for a narrow 
ringlet, and that all the eccentric features in the Laplace Gap also have similar pattern 
speeds (4.72° /day), so there are only seven distinct pattern speeds in Table 01 one for each 
gap/ringlet. The average difference in pattern speed between adjacent features is 0.06° /day, 
with a standard deviation of about 0.01°/day, which is comparable to the uncertainty in the 
precession rates of the patterns (outliers being the 0.04° /day difference between the pattern 
speeds for the Herschel and Russell Gap inner edges and the 0.08° /day difference between 
the Jeffreys and Kuiper Gap inner edges). Also, with the exception of the almost-circular 
Kuiper Gap inner edge, there is an almost monotonic decrease in the amplitude of these 
patterns from 27 km at the Huygens Ringlet to 2 km at the Bessel Gap. This regularity 
hints that these eccentric edges may be controlled by a series of closely spaced resonances. 
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Fig. 6. — Radial excursions of the inner edge of the Barnard Gap, plotted as a function 
of five times the longitude relative to Prometheus. Note that the data are reasonably well 
organized in this longitude system, consistent with the edge being influenced by the 5:4 
resonance with Prometheus. 



This idea will be explored in more detail below after a consideration of the structure of two 
resonantly-controlled edges. 



3.3. Barnard Gap Inner Edge 

The Barnard gap inner edge is a special case because it is the only inner edge of a 
gap in the Cassini Division besides the B-ring edge that cannot be fit to a simple eccentric 
model. All the other non-circular, non-eccentric edges are either on ringlets within the gaps 
(Herschel and Laplace) or at the outer edges of gaps containing such ringlets (Huygens 
and Herschel). Furthermore, the mean radius of the Barnard Gap's inner edge is 120303.7 
km, which is very close to the predicted location of the Prometheus 5:4 Inner Lindblad 
Resonance at 120304.1 km. Thus it is reasonable to expect that the shape of this edge is 
described by the following function: 

r = r D - A* COs(5(A - Xprometheus))- (4) 

Figure [6] plots the radial excursions of this edge as a function of 5(A — Xp r0 metheus) , and 
indeed shows that most of the data can be roughly described by the above functional form. 
There are some clear deviations from the expected pattern, most noticeably some points 
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with large positive excursions near 5(A — Xp r ometheus) = 20°, and a possible phase shift 
in the data relative to the model. Similar deviations are seen in the B ring edge as well 
(see below), and may reflect complications in the dynamics of resonantly controlled edges. 
In spite of this, it does appear that the Barnard Gap is strongly influenced, and probably 
controlled by, the 5:4 Lindblad resonance with Prometheus. 



3.4. B-ring Outer Edge Observations 



The outer edge of the B ring h as long been known to be s trongly affected b y the 
Mimas 2:1 Inner Lindblad Resonance (jGoldreich &; Tremainelll978l : ISmith et al.lll982j ). The 
radial excursions of a particle's orbit near this resonance can be described by the following 
equation (see Section 4 below): 



r = r a - A * cos (2 (A - X Mimas )) (5) 

Analyses of the Voyager data showed that the above expression provided a reasonably 
accurate descriptio n of the outer edge of the B ring in 1981 with an amplitude A ~ 75 km 



(jPorco et al.lll984l ). However, as shown in the top panel of Figure the VIMS occultation 
data do not fit this simple model so clearly. The measurements obtained in 2008 do possess 
a clear m = 2 structure with an amplitude of 75 km, but there is a significant phase offset, 
such that the minimum radius does not occur where A ~ X-Mimas but instead lags behind 
this point by about 40° in phase (i.e. one of the two minima falls ~20° in longitude behind 
Mimas). Furthermore, the data taken prior to 2008 fall well away from the curve described 
by the 2008 data. The radial differences between the 2005-2007 and 2008 data are as much 
as 100 km, which is far too large to be explained by pointing errors, time variability in 
the pattern's amplitude or some additional (smaller amplitude) perturbations in the edge 
position. Instead, the differences between the 2005-2007 data and the 2008 data are best 
explained by a time variation in the orientation of the pattern relative to Mimas. Indeed, 
if we assume the pattern moves at a speed ~ 0.045°/day slower than Mimas' mean motion, 
then the 2005-2007 data are well aligned with the 2008 data (see bottom panel in Figure [7]). 
It therefore appears that between 2005 and 2008 the phase of this m = 2-symmetric pattern 
has drifted backwards relative to Mimas 

This trend in the pattern's orientation over the last few years can also be detected 
by comparing the VIMS occultation data obtained in 2008 with various measurements of 
the B-ring edge made by other instruments earlier in the Cassini Mission. Data from 12 
occultations of the radio signal from the Cassini spacecraft (French et al. in prep) confirm 



1 The long-term average mean motion of Mimas is 381.994509 ±0.000005° /day |Harper fc Tavlorll 19931 ) ■ 



but the 4:2 resonance with Tethys results in a slow variation in the mean motion with an amplitude of 
0.011° /day and a period of 70.8 years. Although this variation is too small to account for the observed in 
the B-ring edge distortion, it does lead to to an instantaneous mean motion at epoch (Day 2005-195) of 
381.984° /day, close to that derived from the JPL SPICE ephemeris. 
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Fig. 7. — Radial excursions of the outer edge of the B ring. The top panel shows the data 
plotted versus the difference in longitude relative to Mimas. The middle panel shows the 
data assuming a constant pattern speed equal to Mimas' mean motion, and the bottom 
panel shows the data organized using a slightly slower pattern speed that provides the best 
fit to the data. Note that data from all 48 occultation cuts in Table [T] are used in these 
plots. 
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that in 2005 the m = 2 radial excursions of the B-ring edge led Mimas by 70 ° in phase 



(or 35° in longitude) . A similar result was obtained by ISpitale Pored (|2006l ) , based on 
Cassini imaging sequences in 2005. The pattern has therefore drifted backwards relative to 
Mimas by ~ 55° in longitude (or 110° in phase) between 2005 and 2008, consistent with 
the drift rate obtained using the VIMS data alone. 

A steady drift in the m = 2 pattern relative to Mimas would be extremely surprising, 
given that this pattern is supposed to be generated by gravitational perturbations from that 
moon. It seems more reasonable that the orientation of the m = 2 pattern instead librates 
relative to the moon on a timescale that is long compared to the Cassini mission to date. 
We therefore posit that the m = 2 structure of the B ring edge can be described by the 
following equation: 

r = r Q - A* cos(2[A - Xuimas - 4>L* sin(Lt - L )]), (6) 

where the last term in the cosine argument describes a slow libration of the longitude of 
the minimum radius relative to Mimas with an amplitude of 4>l and a period of 2tt/L (We 
shall see below that the amplitude A also appears to be time-dependent). 

Estimating (f>L and L from the Cassini VIMS data alone is difficult because the libration 
period appears to be significantly longer than the observation arc. Earlier measurements 
of the B-ring edge orientation back to the Voyager missions could be useful, but these data 
are rather sparse, so before we consider using those data, it is useful to first place some 
constraints on the libration parameters using the Cassini VIMS data alone. 

Since the pattern probably aligned with Mimas sometime during 2007, we can make the 
crude approximation that during the entire Cassini mission to date sin(Lt — 9l) — (Lt — 0L), 
in which case: 

r ~ r Q - A * COs(2[A - (n M imas + <\>iX)t + <AzA - ^Mimas(t = 0)]), (7) 

where nMimas — 381. 99° /day is the current mean motion of Mimas . In this case, the same 
basic procedure described above for finding the pattern speeds of the eccentric edges can 
be used to estimate the value of nMimas + 4>lL that best fits the data. Using all the VIMS 
occultation data the best fit pattern speed for the m = 2 component of the B-ring edge 
is found to be 381.945° /day, or about 0.045°/day slower than Mimas' mean motion (see 
Figure [7] bottom panel; using only the quality data or including the radio science data does 
not change this result much). Thus the product \<PlL\ must be 0.045°/day (where <pL is 
measured in radians and L is measured in degrees/day). 

To turn this estimate of \4>lL\ into a constraint on L, note that (pi cannot have any 
value. First, (f>L cannot exceed 90° or else the pattern would be circulating instead of 
librating, which seems unlikely. Also, (f>L is probably at least 35° or 0.6 rad, since the 
pattern is offset by that much in 2005. These considerations suggest that L most likely 
lies somewhere in the range between 0.03°/day and 0.10°/day, corresponding to a libration 
period of between 10 and 30 years. 

Additional rough constraints on (f>£ and L can be obtained by recalling that analyses 
of the Voyager data showed that the m = 2 pattern was aligned with Mimas in 1980-1981 
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(jPorco et al.lll984l ). This suggests that zero-crossings of the libration occurred in 1980-1981 
and 2007, which would imply the libration period is a submultiple of 52 years, or that L 
is an integral multiple of 0.02°/day. The Cassini-based estimate of \4>lL\ and the Voyager 
data can be satisfied if L ~ 0.04°/day, 0.06°/day, 0.08°/day, etc. (with 4> L ~ 60°, 40° and 
30°, etc. respectively). 

Guided by these constraints, we performed a more comprehensive investigation of the 
B ring edge using a combination of occultation data obtained over the last 30 years. This 
combined data set includes: (1) the 48 VIMS occultation cuts listed in Table Q] (2) 12 
occultati on cuts from 2005 obtained by the Cassini radio-science experiment (RSS) re- 
ported inlFrench et al.l (|2010l ) and kindly provided to us, ( 3) the Voyager 1 RSS occultation 



(Tvle: 



1983 



et al 



1989 
1991 ( 



198 



AM 



19831 ). (4) the Voyager 2 PPS occultation (jLane et al 



1982 



Esposito et al 



mu ltiple data from the ground-based 28 Sagittarius occultations from July 



French et a 



Elliot et al 



1993), and (6) occultations obse rved by the Hubble Space Telescope in 



re- analyzed by lFrench et al.1 (|2010l ). 



1993) and 1995 (IBosh et al.ll2002l ). All of the pre-Cassini data have been 



To determine which combination of libration parameters could best fit these data, 
we computed the Pearson's correlation coefficient p between the radial excursions and the 
parameter cos (2 [A — \Mimas — 4>L * sin(Lt — 9l)]) for different values of (f>z, L and 9l- 
This correlation coefficient will be maximal with the model parameters that best match the 
actual motion of the edge, so it provides a convenient way to search the parameter space 
for likely solutions. This method also has the advantage that it is relatively insensitive to 
any time variability in the amplitude A of the pattern (see below). 

Based on this analysis, we found that solutions with L < 0.05° /day were strongly dis- 
favored. This is primarily because all of the pre-Cassini measurements show negative radial 
excursions, and there was no way to align both the 28 Sgr and Voyager measurements with 
the minima in the Cassini data unless L > 0.05° /day. Furthermore, we found that fre- 
quencies near 0.058°/day and 0.117°/day better organized the data than other frequencies. 
Note that both these frequencies are part of the acceptable series of values that the above 
simplistic analysis suggested would be compatible with both the Voyager and Cassini data. 

Unfortunately, the available data do not clearly favor one of these two solutions for the 
B-ring edge over the other. The top two panels in Figure [8] show the two best-fit solutions, 
with libration frequencies 0.058°/day and 0.117°/day and libration amplitudes of (pL = 52° 
and 32° (in both cases the most recent zero-crossing time was in mid-2007). The scatter 
of the data about the mean curves is not significantly different in these two cases, so both 
solutions are equally good in this respect. Worse, the earlier data are unable to rule out 
the possibility that the m = 2 pattern could be drifting backwards at a constant speed 
instead of librating. The bottom panel of Figure [8] shows the data plotted using the best-fit 
constant pattern speed, which turns out to be 381.932° /day, or about 0.05°/day slower than 
Mimas' current mean motion. The scatter in the data for this solution is not much larger 
than it is for the librating solutions. 



Future observations may eventually provide a way to discriminate between these pos- 
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Fig. 8. — Radial excursions of the B-ring edge, including data from Cassini-VIMS and 
RSS experiments, Voyager observations and various earth-based occultation observations. 
The top two panels show the data as organized using a librating longitude system with 
libration frequencies of 0.058°/day and 0.117°/day, respectively. The bottom panel shows 
the data organized using a steady drift rate of ~ 0.05°/day slower than Mimas' mean motion 
(as before, times are measured relative to an epoch of 2005-195T02: 12: 13.557 or a Cassini 
spacecraft clock time of 150000000). The available data are clearly insufficient to distinguish 
between these different models, given the very large residual scatter 
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Fig. 9. — Radial excursions of the B-ring edge measured by the VIMS and RSS experiments 
in 2005 (top) and by VIMS in 2008 (bottom). In both cases, the data can be reasonably 
well fit by simple sine curves. 



sible solutions. However, some of the difficulty in determining the correct model for the 
motion of the m = 2 pattern stems from the comparatively large scatter in the measure- 
ments with respect to any of the above solutions. In all cases, the average amplitude of the 
pattern is around 60 km, but rms residuals are roughly 20 km, which is much larger than 
the ~ 1 km measurement errors in each of these data sets. This suggests that the B-ring 
edge is not just a fixed-amplitude m = 2 pattern that librates relative to Mimas, but has a 
more complex shape, perhaps with additional perturbation modes. 

To explore this possibility, let us take a closer look at two particular subsets of the data: 
the VIMS data from 2008 and the combined VIMS and RSS data from 2005. Each of these 
subsets consists of a reasonably large number of occultation cuts (31 and 16, respectively) 
that were not only taken over a sufficiently broad range of longitudes relative to Mimas that 
we can estimate the shape of the m = 2 pattern, but also obtained in a sufficiently short 
period of time that we can ignore the libration of the pattern with respect to Mimas. 

The best-fit m = 2 pattern for the 2008 data has an amplitude of 75.3 km and a phase 
of -39.4° relative to Mimas (i.e. the pattern lags Mimas by 19.7° in longitude), while the 
best-fit pattern for the 2005 data has an amplitude of 48.7 km and a phase of +70.0°. Thus 
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Fig. 10. — Residual radial variations of the outer edge of the B ring, after removing the 
m = 2 mode pattern associated with the Mimas 2:1 resonance. The top and middle panels 
show the 2008 and 2005 residuals, respectively, each organized using the best-fit pattern 
speed. Note the pattern speed for the 2008 data is close to the apsidal precession rate at this 
location, while the 2005 data is best fit by a 20% slower pattern speed. The bottom panel 
shows the 2005 data plotted using the same pattern speed as the 2008 data, to illustrate 
that the data are relatively well organized in this case as well. As before, all times are 
measured relative to an epoch of 2005-195T02:12:13.557 or a Cassini spacecraft clock time 
of 150000000. 
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the amplitude of the m = 2 pattern varies with time as well as the phase (see Figure[9j). This 
is not entirely unexpected, given that a libration can often be modeled as a combination of 
free and forced terms, which naturally leads to coupled variations in both the phase and the 
amplitude of the pattern (see Section 4 below). Indeed, some of the observed scatter in the 
data in Figure [8] is likely due to unmodeled changes in the amplitude of the m = 2 pattern. 
Even so, the variations in the amplitude of the pattern cannot be the only source of scatter 
in these data, as each of the 2008 and 2005 data sets alone show significant scatter around 
the mean m = 2 pattern, with individual residuals up to ~40 km (see Figured]). 



A previous analysis of Cassini images by ISpitale fc Pored (|2006l ) suggested that addi- 



tional perturbations, possibly m = 3, might influence the B-ring edge's shape. Motivated 
by this result, we took the 2008 and 2005 data sets, removed the best-fitting m = 2 pattern 
from each of them, and fitted the residuals to m = 0,1,3,4,5, and 6 patterns. For each 
non-zero value of m, we ex amined a range of p attern speeds in the vicinity of the expected 



speed of a "normal mode" ([French et al.lll99ll ) with this m— number: 



rip = [(m — l)n + w]/m. (8) 
For m = 0, the pattern does not rotate, but the whole ring osc illates radially at th e epicyclic 



frequency k = n — w, as observed for the Uranian 7 ring ([French et al.lll986l ). As with 



eccentric features in the Cassini Division, we computed the amplitudes: 

1 n 

- 5ri * cos[m((l - 5 m0 )Xi - Q p 5ti)} (9) 



atR 

n . 
i=i 



aj 

n 
i=l 



1 n 

E5ri * sin[m((l - 5 m0 )\i - tt p 5ti)} (10) 



where Sri are the residual radial variations in the edge positions after removing the m = 2 
pattern. As before, if there is a pattern in these radius measurements with a given m and 



Op, then w a R + aj should have a maximum at the appropriate values for those parameters. 

For both the 2005 and 2008 data sets, the clearest maxima were obtained with m = 1. 
Indeed, as shown in Figure \W\ the residuals in the 2008 data are reasonably well fit by 
a precessing Keplerian ellipse with an amplitude of about 18 km and a pattern speed of 
5.08°/day (with an uncertainty of ~ 0.05°/day determined by the ~6-month observing arc). 
This pattern speed is very similar to the expected apsidal precession rate of 5.061°/day for 
a Keplerian orbit at this radius. For the 2005 data, the best-fit solution has a slightly larger 
amplitude of about 23 km and the best-fit pattern speed of 4. 19° /day, 17% slower than the 
predicted value of w. However, given the short time-baseline covered by the 2005 data and 
the limited number of data points available, a pattern speed of 5.08°/day (consistent with 
the 2008 observations) can also organize the 2005 data reasonably well (see bottom panel 
of Figure [TO]). An attempt to fit a single, coherent m = 1 perturbation to both data sets 
yields multiple equal-quality fits at an array of pattern speeds separated by 360° /3 years 
~ 0.3°/day. The available data are therefore too sparse to permit us to derive a more precise 
model of this perturbation or an accurate measurement of the pattern speed. However, the 
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existence of similar m = 1 patterns in both the 2005 and 2008 data sets strongly suggests 
that this is a "permanent" feature of the B-ring edge. 

The radial variations in the location of the B-ring edge therefore have at least two com- 
ponents: an m = 2 perturbation forced by the strong Mimas 2:1 Inner Lindblad Resonance, 
with an average radial amplitude of ~60 km and an orientation which librates (or circulates) 
with respect to the direction to Mimas, and an m = 1 Keplerian ellipse with an amplitude 
of ~ 20 km which freely precesses under the influence of Saturn's oblate figure at roughly 
5°/day. Even this rather complicated model provides a relatively poor fit, in comparison 
with those for the eccentric edges in the Cassini Division, with rms residuals of 10-20 km. 
We suspect that this is due to deficiencies in our model of the m = 2 perturbation, whose 
libration is still not accurately modeled (e.g. we have not yet implemented a model for 
the amplitude variations), but it is also possible that additional perturbations are present. 
Efforts to derive more detailed and accurate models of the B-ring edge are underway and 
will be presented in future work, but already the preliminary results presented here hint at 
a close connection between the B-ring edge and the Cassini-Division gaps. 



4. Theoretical expectations for the B-ring outer edge 

Before delving into the possible connections between the B-ring edge and the Cassini- 
Division gaps, it is useful to first take a closer look at the complex behavior of the B-ring's 
outer edge. In particular, the above data allow several different solutions for the motion of 
the m = 2 component on this edge, all of which have the m = 2 pattern drift or move relative 
to Mimas. We would like to establish whether any of these solutions are plausible in terms of 
the local dynamica l envir onme nt. Since a detaile d model along the lines of those developed 



m iBorderies et al.1 (119821 ) and lHahn et al.l (|2009l ) is beyond the scope of this paper, we will 
instead examine the behavior of isolated ring particle orbits in the vicinity of the Mimas 2:1 
inner Lindblad resonance (ILR). While not conclusive, these simpler calculations do at least 
demonstrate that the observed variations in the amplitude and orientation of the m = 2 
pattern are not wildly inconsistent with theoretical expectations. 

The following discussion is based on that in Chapter 8 of Murray and Dermott (1999), 
but couched in terms of the physical coordinates, r and A and the standard orbital elements, 
a, e and w, rather than a Hamiltonian formalism. 

At a first-order ILR, the resonant argument is given by 

fiLR = (m - 1)A + zd - m\ s (11) 

where m is a positive integer, A and A s refer to the mean longitudes of the test particle and 
satellite, respectively, and zu is the longitude of pericenter of the test particle's orbit. At 
the so-called exact resonance, 

dt = {m - l)n + OT SCC - mn s = (12) 
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where n is the orbital mean motion and w sec is the apsidal precession rate of the test particle 
due to non-resonant perturbations (chiefly the planet's oblateness). For a 2:1 ILR, m = 2 

and 7i r es = 2n s — w sec . 

An ensemble of test particles which share common values of semimajor axis a, eccen- 
tricity e and <^ilr, but different instantaneous values of A and w, will define a streamline 
given by 

r(A, t) « a [1 — e cos(A — zu)] , (13) 

or in terms of 9?ilr: 

r(A, t) w a [1 - e cos (m(A - A s ) - (^ilr)] • (14) 
This is the same as Eq.([6]) above, with A = ae, r' Q = a and 

piLR = 2fosin(Lt-0 L ). (15) 



Following Murray and Dermott (1999, see their Eq. (8.26) and Table 8.5), the corre- 
sponding time-averaged disturbing function of the satellite (to lowest order in eccentricity 
and for zero inclination) is given by 



7^-ilr 



Gm, 



/(a)ecos <£ilr 



(16) 



where m s is the mass of the satellite, a = a/a s and for m = 2 the function f(a) = —0.75/a 
comes from evaluating the disturbing function. The Lagrange equations which describe the 
resulting perturbations in the test particle's orbital elements are: 



dn 
~dt 



dw 
~dT 



■ -3(m - l)/3n 2 esin<^iLR, 

de a ■ 

— = -pnsmtpiLR, 

at 

-- -fine' 1 cos ipi LR + w sec , 
-j t = --pne cos ipj LR , 



(17) 

(18) 
(19) 
(20) 



where we have introduced the dimensionless resonance strength (5 = —(m s /Mg)af(a), with 
Ms the mass of the planet. The quantity e is the longitude at epoch, defined by the usual 
expression A = e + nt. These are equivalent to Eqns. (8.28-8.31) of Murray and Dermott 
(1999), except that our definition of the resonant argument in Eq. (|11|) is opposite in sign 
to theirs. Combining Eqns. <\17\ [T9l and I2T)|) . we can derive an expression for the rate of 
change of c^ilr 



diflLR 

dt 



(m — l)n — mn s + w sec 



fin (m — l)(3ne 
~ + 2 



cos ipi LR . 



(21) 



For small values of e, we can safely neglect the last term in brackets (which arises from 
Eq. p}). 
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From Eqns. (|17l and[T8 |) . we see that n and e must vary in phase, with 

dn 

— = 3(m - l)ne, (22) 
ae 

so that we may write, to the lowest order in e 

3 

n ~ n + -(m - l)n e 2 , (23) 

where no is a constant. Substituting this expression for n into Eq. (|2ip . and introducing 
the constant parameter v = (m — 1)uq — mn s + vu sec , we have our final equation for the 
resonant variable, correct to the lowest order in e: 

— — = i/ + -(m - 1) n e cosc^flr. (24) 

ar 2 e 

The frequency parameter i/ is best thought of as a measure of the distance from exact 
resonance, 

v ~ -|( m ~ 1 ) n res f - — — ) (25) 



2 V a 



'res 



where a res and the semimajor axis and unperturbed mean motion at exact resonance, 

defined by Eq. (fT2|) . Note that v > for orbits interior to a r 



'res • 



Equilibrium solutions to Eqns. (fTHI and [24"J) occur where e = eo and <^ilr = V^Oj where 
eo and t^o are constants in time. Such solutions only exist for (fo = or ipo = ir, and where 
the eccentricity satisfies the cubic equation: 

P n 3. 2 2 

i/ = cos<^o - -( m - 1) n oe - (26) 

eo 2 

For positive values of v (i.e. interior to a res ) only a single solution exists, with </?q = 0. But 
for negative values of v (i.e. exterior to a rcs ) up to three branches of solutions exist, labelled 
A, B and C by Murray and Dermott (1999, see their Fig. 8.9). Solution A, with ip$ = 0, is 
an extension of the equilibrium solution for v > 0. Solutions B and C have ipo = tt, with 
< &c < ca- Solution C is an unstable equilibrium, and thus of little physical interest, 
but solutions A and B are both stable. Branches B and C merge at a bifurcation point, 
where dv/deo = 0, 

^ = -^n [3(m-l) 2 /5 2 ] 1/3 , (27) 



and 



3(m - l) 2 



1/3 

(28) 



and cease to exist for larger values of v. 

For e << e c , the equilibrium, or forced eccentricity is given approximately by 

60 ~ P-^ COSL Po, (29) 
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with ipo = for v > (solution A) and ifQ = it for v < (solution B). For e >> e c , ipo = 
(solution A) or 7r (solution C) and 

1/2 



eo 



2M 



3no(w, — l) 2 



(30) 



At v = 0, e = 2 1 / 3 e 



In the context of planetary rings, it is usually assumed that we are in the 'small-e' 
regime, with e « e c and a resonantly- forced eccentricity and phase given by Eq. (|29|) (cf. 
Murray and Dermott 1999, Eq. (10.21) and Fig. 10.9). 

If we now substitute suitable numerical values for the Mimas 2:1 ILR, we find that 
a res = 117, 553.65 km, w^.r. = 5.0613° /day, n res = 2n s — w sec = 758.9277° /day, a pa 0.6336, 



m s /M = 6.60 * 10~ 8 (jjacobson et al.1 120061 ). and j3 = 4.95 * 10~ 8 . The corresponding 
eccentricity and frequency at the bifurcation point are e c = 0.00255 and v c = —0.0221° /day. 
From Eq. (|25[) we have that 



v pa -0.0097° /day 

so that the bifurcation point occurs at a — a res : 
299 km. 



a — a r 



1 km 

+2.3 km with a forced amplitude, ae c 



(31) 



Since the maximum radial amplitude observed at the edge of the B ring, A ~ 75 km, 
or one-quarter of the critical value, we may conclude that the streamlines are indeed in the 
'small-e' regime. We can then use Eq. (|29p to estimate the effective value of v for streamlines 
at the edge of the B ring. If we adopt an average forced eccentricity of (0.50 ± 0.12) * 10~ 3 
(i.e. A = 60 ± 15 km), and note that tp appears to librate around rather than ir, then we 
find 

0n o 



eo 



• cos ip Q = +0.075 ± 0.019° /day 



(32) 



with an effective distance from exact resonance of a — a Tes = —8 ± 2 km. This may be 
compared with the observed mean location of the edge of 117,565.4 km, or a—a res = +12 km. 
The B ring edge thus behaves as if it were located 5 — 10 km interior to the 2:1 resonance, 
whereas it is actually located ~ 12 km exterior to the resonance. A similar conclusion was 
reached bv lPorco et al.l (|1984l ) in their study of this edge using Voyager observations, though 
their value of a — a Tes was less accurate due to uncertainties at that time in the absolute 
radius scale of the rings. 

However, the classical 'resonance width' w res — the range in semimajor axis over 
which test particle orbits on opposite sides of exact resonance will overlap due to their 180° 
difference in ipo — is given by substituting Equation [25] into Equation [29] and solving for 
/2 (cf. Murray and Dermott 1999, Eq. (10.23)) 

nl/2 



Wy 



2a r 



43 km. (33) 

3{m — 1) 

(Note that this is actually 1/2 of Murray and Dermott's expression, which includes the 
radial excursions due to the forced eccentricities.) Within ~ 20 km of exact resonance, 
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therefore, we cannot expect the above test particle model to give realistic estimates of eo 
as a function of a. In an actual planetary ring, the situation will be further complicated 
by gravitational and collisional interactions between the ring particles, which will act to 
prevent the streamline crossing predicted by our simple test-particle model within ±w res /2 
of a rP s- Such a model has recently been develope d for application to th e A and B ring edges 



by lHahn et al.1 (|2009l ) , based on earlier work by iBorderies et al.l ()1982l ) . 



Although dissipative collisions within a real ring might be expected to rapidly damp any 
motion relative to the equilibrium solutions discussed above, our observations of the B ring 
edge strongly suggest either that the resonance angle ¥>ILR is librating about an equilibrium 
value close to 0, or that it is circulating in a retrograde direction (i.e. (dipi^/dt) < 0). The 
equations of motion, (fT71 - l20"|) admit of both finite-amplitude oscillations in 9?ilr and of 
circulating solutions. For a series of phase portraits of such librations the interested reader 
is referred to Fig. 8.10 in Murray and Dermott (1999). (Their dimensionless parameter 5 is 
equal to 3z//f c , while their amplitude parameter $ = ^(e/e c ) 2 .) 

In the small-e limit, (i.e. e « e c ) Eq. (f24"|) reduces to 

d<piLR Pn Q 

—7— ~ v cos Vilr, (34) 

at e 

while the variation in e is given by Eq. (|18p . where we may set n ~ uq. In this limit, the 
variations in <^ilr are dominated by the resonant effects on zu. These coupled equations 
are most readily solved by introducing the conjugate variables 

h = e cos ipi LR , k = e sin <£ ILR , (35) 

in terms of which Eqns. (I18p and (I34D become 

dh , % 

* = -"«=■ «> 

d —= v }i- /3n . (37) 

The solution follows trivially: 

/?rao ^ ( 38 ) 

v 

k(t) =e/sin(i4-0). (39) 

Recall that (3n/\v\ is the forced eccentricity, eo from Eq. (I29h above. The free eccentricity, 
ef and phase angle 6 are arbitrary constants of the motion, set by the initial conditions. In 
the h, k plane, the motion is in a circle of radius ej about the fixed point, (±eo, 0). Interior 
to the resonance [y > 0, equilibrium branch A) the motion is counterclockwise about (+eo, 
0), while exterior to the resonance {y < 0, branch B) the motion is clockwise about (— eo, 
0). In both cases, the angular frequency is simply equal to given by Eq. ([25]) above. 

If ef < eo, the solution describes a librational motion of ip about either (for u > 0) 
or 7r (for v < 0). The maximum libration amplitude is given by tp mSLX = sm_1 ( e // e o)> while 
the instantaneous eccentricity is equal to 

e(t) = Vh 2 + k 2 = [e 2 Q + e) + 2e e/ cos(z/t - 6)] 1/2 . (40) 
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For ef > eo, the angle </?ilr circulates continuously through 2n radians while e oscillates 
between ef — eo and ef + eo- Example trajectories for both cases are shown in Fig. 8.11 of 
Murray and Dermott (1999). 

Although the form of the small-e librations motivated our choice of the model used 
to fit the B ring edge above, we note that even this simple test-particle model implies (i) 
that the amplitude of the edge, A = ae should oscillate as the edge librates or circulates, 
and (ii) that the variation in the resonant angle, v?ilr will not be sinusoidal or linear, 
unless ef « eo or e/ >> e Q . The substantial variations observed in A and v?ilr over the 
course of the Cassini mission to date suggest that the latter is unlikely to be true. A more 
sophisticated model might therefore be based on the small-e solution above. 

Despite these limitations, we note that the best-fitting libration frequencies of 0.058° 
and 0.1 17° /day found above are generally consistent with our theoretical estimate of v = 
0.075 ± 0.019° /day in Eq. (132[) . based on equating the observed amplitude to the forced 
eccentricity. 

If, on the other hand, the B ring edge is circulating, then the best-fitting drift rate of 
(j) L L = -0.05° /day (cf. Fig. 8) implies that d<pi LR /dt = 2cf> L L « -0.10° /day, from Eq. (fT5j) . 
or a — a res = +10 km, which agrees fairly well with the observed mean radius of the edge. 
Thus both librating and circulating small-e models seem to be quantitatively consistent 
with our observations of the B ring edge, despite their obvious limitations noted above. 

For completeness, we did examine the large-e case, but these calculations indicate that 
a large-e libration is less compatible with the observed motion of the B ring edge. Not 
only does the predicted forced amplitude, A = aeo exceed 750 km only 5 km outside the 
resonant radius — ten times the maximum observed value — , but the libration frequency 
is quite slow, reaching only 0.024°/day at a — a res = 5 km and 0.033°/day at 20 km. Such 
low libration frequencies are strongly disfavored by the available data. 

5. A possible explanation for the location of the Cassini Division gaps 

The above analysis indicates that the structure of the Cassini Division may be more 
regular than it appears at first glance. Certainly, the Barnard Gap is no longer a mystery, as 
it seems to be held open by the 5:4 mean motion resonance with Prometheus. Furthermore, 
the remaining gaps seem to form a pattern. The Herschel, Russell, Jeffreys, Kuiper, Laplace 
and Bessel gaps all seem to have an eccentric inner edge and most of these also have a circular 
outer edge. (The Huygens, Herschel and Laplace gaps do not have perfectly circular outer 
edges, but these are also the gaps which contain dense ringlets, which may complicate the 
situation.) Furthermore, the 6 eccentric gap edges seem to define a series of evenly spaced 
pattern speeds with a characteristic spacing of about 0.06° /day. Even the Huygens ringlet 
and the m = 1 component of the B ring edge seem to fall in this sequence. Finally, there 
is the suggestive coincidence that the m = 2 pattern on the B ring edge could be librating 
with an angular frequency not too dissimilar from 0.06° /day (although we must caution 
that other solutions are possible). 
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Based on these findings, we have developed a novel explanation for the location of the 
gaps in the Cassini Division. Just as the inner edges of the Huygens and Barnard Gaps 
are established by resonances with saturnian satellites (Mimas 2:1 and Prometheus 5:4, 
respectively), the inner edges of the other gaps are maintained by resonances which involve 
the edge of the B ring. In the following sections, we demonstrate that the gravitational 
perturbations from the radial excursions of the B-ring edge, together with perturbations 
from Mimas, can give rise to terms in the equations of motion that would support the 
formation of eccentric edges at their observed locations in the Cassini Division. 

Note that for the purposes of this theoretical study, we will assume that the orientation 
of the m = 2 pattern on the B-ring edge librates with a frequency of ~ 0.06° /day, which 
matches the typical spacing between the pattern speeds of eccentric features in the Cassini 
Division. This particular solution is clearly the one where the B-ring-edge patterns would be 
most likely to generate a series of structures like the Cassini Division gaps. Thus an analysis 
of this case will establish whether such a model has any hope of working. Of course, our 
current understanding of the data admits the possibility of different libration frequencies 
or even circulation of the m = 2 pattern. If one of these alternate solutions turns out 
to be correct, then the connection between the B-ring edge and the Cassini Division may 
still exist but be more complicated, and the simpler analysis presented here could still help 
clarify this relationship. 



5.1. Qualitative Frequency Studies 

To determine if any interesting resonances in the Cassini Division could be generated 
by the B-ring edge, one must write down the gravitational potential associated with the 
observed structure of the B-ring edge, compute the appropriate terms in the disturbing 
function, and determine if they could drive structures like those seen in the various gap 
edges. In this case, however, the data described above already provide hints of what terms 
in the disturbing function could be involved in generating the observed edges. Therefore, 
prior to exploring the dynamical environment of the Cassini Division in detail, we will 
first take a more qualitative look at the situation in order to clarify which terms in the 
gravitational potentials are likely to be relevant. 

The m = 1 patterns on the relevant Cassini Division edges move at speeds close to the 
expected local apsidal precession rates given Saturn's oblateness (See Tabled]). Thus each 
edge can be thought of as a collection of particles on freely-precessing eccentric orbits whose 
pericenters are all aligned to produce a coherent structure. Assuming the orientation of the 
m = 2 pattern in the B-ring edge librates at ~ 0.06° /day, this alignment in pericenters 
would appear to occur at places where the apsidal precession rate w has the following 
values: 

W = TD B - jL (41) 

where wb is the apsidal precession rate of the B-ring edge (also the pattern speed of the 
m = l component on that edge), j is an integer and L is the libration frequency of the 
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m = 2 pattern in the B-ring edge. Since zd as a function of a is determined primarily by 
the higher-order components in Saturn's gravitational field, Equation [41] implicitly specifies 
the locations of a series of regularly-spaced resonances. At these locations, the following 
quantity is approximately constant: 

<p = zd — zdb + jLt (42) 

We propose that there is a term in the equations of motion that tries to maintain ip near 
some value ip a . Such a term will act to align orbital pericenters and produce a coherent 
structure on each edge that is stable against small perturbations. In order for this to work, 
the equation of motion of this resonant argument needs to have a term of the form: 

^~-/ o 2 sim>-^ ). (43) 

In this particular situation, the above expression can be re-written as (assuming zdb and L 
are constant): 

dP ZD 

-p- ~ -fl * sm(w -zd b + jLt - ip ) (44) 

Provided such a term exists in the equations of motion, then the pericenters of individual 
particle orbits could become aligned, forming a coherent m = 1 pattern that moves around 
the planet at the local precession rate established by Saturn's oblateness, consistent with 
the observations. 

Lagrange's Planetary Equations relate time derivatives in orbital elements like zd to 
derivatives of the disturbing function (e.g. dw/dt cx dTl/de). Thus, in order to determine 
if if = zd — zdb + jLt is a proper resonant argument that can produce a coherent structure, 
we need to find a term or a combination of terms in the disturbing function that are 
proportional to s'm(w — zdb + jLt). 

The question now is whether such terms are likely to appear in the disturbing function 
generated by the B-ring edge. We can model the B ring edge as a collection of mass ribbons, 
each of which has the characteristic shape 

r = r Q - di cos(0 - w B ) - d~2 cos(2[0 - Xuimas - 4>l\) (45) 

where <f> is the azimuthal angle. <Pl = 4>L * sin(Lt — 6l) is the time-variable azimuthal 
offset due to the edge's libration (cf. Equation [6] above). Note the term c?2 may also have a 
time- variable component (c^ = (fa[l ±5^ cos(Li — Ql)]) because the amplitude of the m = 2 
pattern is apt to change over the course of the libration cycle. In fact, as discussed in 
Section 4 above, the librations in the amplitude and the phase are likely to be coupled in 
such a way that neither oscillation is purely sinusoidal, but for the purposes of this analysis 
we will ignore such complications. 

Such mass ribbons will produce a relatively complicated potential, but in general we 
expect terms in the disturbing function to contain second-order terms of the following form 

Tiki OC 4 fc '4 Z ' COs(fc(A - ZD B ) + Z[2A — 2\ M im,as - 2</>l]), (46) 
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where k and I are various integers and A is the perturbed particle's mean longitude (assuming 
the particle's orbit is nearly circular). We can quickly see that any term of this form that 
contains 4>l will also contain \Mimas, so no term of this form will have the form required to 
produce something like Equation 1441 This suggests that the B ring edge alone will be unable 
to produce the resonances required to explain the eccentric edges in the Cassini Division. 

All is not lost, however, because there is a possibility that the perturbations from the 
B ring edge, acting together with perturbations from the nearby and very powerful Mimas 
2:1 Inner Lindblad Resonance, can generate the desired terms in the disturbing function 
through a sort of three-body resonance. In such a situation, the time-derivatives of orbital 
elements can be proportional to the products of disturbing functions, which means we will 
get products of the IZti above with the term responsible for the Mimas 2:1 resonance: 

H-Mimas OC m M COs(A + W - 2X M imas), (47) 

where rriM is the mass of Mimas. This yields terms of the following form: 

TZ' k i oc m M df^2 sin((A + w - 2\ M imas) + k(X - w B ) + l[2\- 2\ M imas - %4>l\) (48) 

When k = 1 and I = — 1 the term has the promising form: 

72-1 _! oc mMd\d2 sin(w — zub + 24>l)-, (49) 

or, equivalently 

TZ'i _i oc raud\d2 sin(ro — wb + 24>i sin(Lt — 6l))- (50) 

This term is proportional to something like sin(A * t + B * sin(C * t)). Note in particular 
that in this case (pi — 50° or nearly 1 radian (cf Fig. |Sj), so the coefficient B is actually 
quite large. Figure [TT1 shows a fourier transform of a function of the form sin(t + sin(.Olt)), 
which illustrates that this sort of function actually has multiple periodic components (the 
amplitudes of these different components can be evaluated using Bessel functions, see Gray 
&i Mathews 1895). The above term in the disturbing function can therefore be expressed 
as the following series: 

Tt'i-i oc Cj sin(ro — wb + jLt — j8jj), (51) 

3 

which is exactly the form which we are seeking (cf. Equation 1441 above) . A single disturbing 
term of the form given in Equation [49] could therefore potentially lead to the formation of 
multiple, evenly-spaced features, like the gaps in the Cassini Division. 

Given all this, it seems reasonable to look for three-body-like resonances in the Cassini 
Division where one of the bodies is Mimas and the other is the B-ring edge. Furthermore, 
we will be particularly interested in terms in the disturbing function of the B-ring edge that 
have the form: 

72-1 _i oc did 2 sin([A - w B ] - [2A - 2X Mimas - 2<f> L ]) (52) 
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Fig. 11. — Fourier transform of the function x = sin(i + sin(.Olt)), showing the multiple 
frequency components with comparable amplitudes. 



5.2. Evaluating the Disturbing Function 

The above analysis indicates that we will need two pieces of the disturbing function: 
(1) The part of the disturbing function due to Mimas that goes like cos(A + w — 2\Mirnas) 
and (2) the part of the disturbing function due to the B ring that goes like cos(A + wb — 

The relevant part of the Mimas disturbing function is given by Equation [16] above: 

T^Mimas = GmM f(a/a M )e cos(X + w - 2\ M im as ) (53) 

where ttim and au are Mimas' mass and orbital semi-major axis, respectively, while a 
and e are the particles' semi-major axis and eccentricity, (recall that for a 2:1 resonance 
f(a) = -0.75/a). 

Of course, extracting the relevant bit of the ring's disturbing function requires more 
effort. Again, consider a mass ribbon whose position is described by the following equation: 

r = r a — 5r = r Q — d\ cos(0 — wb) — d 2 cos(2</> — 4>lm) (54) 

where for simplicity of notation we have introduced the term: 4>lm = ^Xuimas + 2(f>L- 

The disturbing function for a small mass element dm of this ring on a particle at radius 
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r' and longitude A is given by: 



^dm = Gdm 



1 OO , i 

1 / t \ 1 , , . r 



cos V> 



(55) 



where ^ is the difference in longitudes between the mass element and the particle of interest 

„„ A „/. i \ 



so 4> = ip + A. 



The total disturbing function for the ring is then this expression integrated over all tp: 



Turing — 



Gm r f 2n 

~2^~ 



Jo L i=i 



cos V> 



(56) 



where m r is the total mass of the ribbon, which for now is assumed to be uniformly dis- 
tributed in (f>. 

The terms we are interested in here 



(Sr 



are proportional to (Sr) 2 = (r — r a ) 2 because: 
2d\d 2 cos(V> + A — wb) cos(2?/> + 2A — <Plm ) + •••• 



(57) 
(58) 

. . . (59) 
So this is the lowest order term in d±, d 2 that contains the desired frequency term. 

Re-writing lZ r ing so that Sr is explicit, we find: 



(Sr) 2 = d\d 2 cos(V> + A + wb - 4>lm) + 
(Sr) 2 = d\d 2 cos(A + w B - 4>lm) cos ip + ■■■■ 



ring 



Gm r 
2irr' 



1 = 1 V ' 



2ir 



(r Q - Sr)' 



■ cos ipdip (60) 



Expanding both terms and keeping only terms proportional to (<5r) 2 ,we obtain: 



Gm 



ring 



2-nr'r 



712 



1(1-1) fr Q \l 



2n 



(Sr) 2 P;(cos ip)dtp 



I i-2-k 



3Gm r r' 

2^ ,„ 



(Sr) 2 cos i/jdip 



(61) 

While both parts of this function are non-zero, the second part is a single term while (as we 
will see below) the first part is a slowly converging series. Thus at this point we will drop 
the second term and only keep the parts of the first term containing the element explicitly 
listed in Equation 1591 

r>2vr 



ring 



Gm r did 2 , t (r Q \ l 
-j-Z- cos(A + w B - 0lm) 2^2\?. 

1=2 



2iTr'r^ 



(I — 1) cos ipPi (cos ip)dip 



(62) 



It turns out that for large values of I the term in the square brackets becomes 4 for 
odd I and for even I (see Fig. I12[) . Since r a /r' ~ 1 in the Cassini division, the series slowly 
converges, so we can actually approximate this term as 2. In this approximation: 



, Gm r did 2 , ~ ^./^V 

n ring ~ g , 2 COs(A + W B ~ <pL M ) 2^ l [zj) > 
° 1=2 



(63) 
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Fig. 12. — The term J 27r (n — 1) cos if)P n (cos if})dip, evaluated for different values of n. The 
bottom branch corresponds to even n, where the integral is zero by symmetry. The top 
branch corresponds to odd n, which asymtotically approaches 4. 



and since lx l = 1/(1 



x) 



Gm r r 



12 



d\d 2 



rm9 27rr' rl (j J - r )'- 



' C0S(A + TUB - 4>LM) 



(64) 



Finally, we assume in this case that the particle is on a nearly circular orbit, so r' ~ a 
so the final expression for this function is: 

Gm r a 2 d\d,2 



■ cos(A + w B - 4>lm) 



(65) 



2ira rl (a — r 

where the superscript r indicates that this perturbation is due to the radial distortions in 
the B ring edge. 

This expression gives the disturbing function for a single ribbon of mass total m r . This 
perturbation from the entire outer part of the B ring can be approximated as a series of 
such ribbons with different v alues of m r , r , <ii and d 2 - Since the structure of the B ring is 
still uncertain at this point ( Hahn et al.ll2009l ). for simplicity we will here assume that d± 



and di decrease linearly with radius towards the planet: 

, r n -Ri 
d\ = Di 



do = D< 



R Q — Ri 

r Q - Rj 
Ro — Ri 



(66) 
(67) 
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where R a is the mean radius of the B ring edge, D\ and L>2 are the values of d\ and di 
at the edge, and Ri = R Q — W is the assumed radius where the amplitudes of the radial 
excursions in the mass ribbons go to zero. Furthermore, let us assume that the mass of 
each ribbon is: 

m r = 2irar dr (68) 
where a is the (unperturbed) surface mass density of the outer B ring. 

If we insert the above expressions for m r , d\ and cfo and integrate over all values of r Q 
between Ri and R Q , assuming that W = R Q — Ri is much less than R Q and a — R Q , we find 
the disturbing function from the entire B ring is given by: 

GaW a Lhth n _^ ~, , 

K Ring ^ (a-R) 2 COS ^ WB ~ ^ LM > ^ ' 

At this point, we should note that the above model of the B ring might be too simplistic, 
in that it assumes the mass ribbons are of constant mass. It is of course possible that the 
mass density of the ribbons p also varies with longitude: 

p=— L [l + Hi cos(0 - w B ) + fJ-2 cos(20 - <Plm)}, (70) 

where (ii and [12 are fractional mass density variations along the ribbon. If particles do 
not collide with one another, then the mass variations p,\ and ^2 could be derived from 
the gradients in the radial variations and orbital parameters. However, in practice the 
outer part of the B ring is very dense, so collisions are very important. Indeed preliminary 
investigations show that the local optical depth of the B ring can vary by over a factor of two. 
While we expect these density variations to have some connection with the radial positions 
of the outer edge, at the present moment we will avoid attempting to estimate p,\ and ^2- 
Regardless of the numerical values of these parameters, such mass anomalies would mean 
that terms of order [i\d2 and would also produce terms with the desired frequency. 
Following the same procedures as described above, one can show the disturbing function for 
a ring with both radial excursions and mass anomalies is (ignoring the possibility of phase 
shifts between the mass anomalies and the radial excursions): 



K ^ * —R- Ja-R^f ~ ^^~R -^^-R ) C ° S(A + WB ~ (71) 



5.3. Numerical Values of the Disturbing Function Terms. 

Having obtained these two terms, it is worth evaluating their coefficients so we can 
make some quantitive comparisons of their strength: 

For T^Mimas this is relatively easy, since: 

f(a/a M )ecos(\ + w - 2\ M imas) (72) 

AM 
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or, equivalently: 

f(a/a M )(ae) cos(A + w - 2X Mimas ) (73) 

Given rau — 4 * 10 19 kg, and (a/ o,m) f (a/ Q-m) — —0.75 for a 2:1 resonance, and assuming 
a ~ 118, 500 km and ae ~ 5 km (appropriate for the Russell Gap), we find: 

_ 4 2 2 / ae \ / 118, 500/cm\ 2 
K-Mimas = -7.1 * 10 m/s (jjT^JI I cos(A + - 2X Mimas ) (74) 

For lZm n g there are more uncertain terms. First considering the simpler case of pure 
radial excursions in the ring (Equation [69]) . 

GaW a Dx£> 2 - , . . 

ft-Ring - 3 ~f{~ (a — R ) 2 COS ^ +WB ~ ^ LM > ( 75 ^ 

Again, assume a ~ 118,500 km for reference. Also, based on our measurements, we will 
assume R Q ~ 117, 500 km , D\ ~ 20 km and D2 — 60 km. For rough purposes, we will use a 
surface mass density of a ~ 100 g/cm 2 and W ~ 100 km, though these parameters are less 
well constrained. Substituting these numbers into the above expression gives the following 
result: 



B r w .. ~ 2.7 * 10 6 m 2 /s 2 * I : I I ] I ^ ] I — ] I — I I I cos(A + H7 B - cj> LM ) (76) 

- R "* 9 \W0g/cm 2 J VlOOfcm/ \ 118, 500fcm / \ 20km J \60km J \{a- R a ) J 

Thus the perturbations from the B ring are about 200 times weaker than those from the 
2:1 Mimas resonance. While this is a rather large factor, it is not many orders of magnitude. 
Also note that the perturbations from the ring could be even larger if there are significant 
variations in the mass of the ribbons with longitude. 



( a \ ( w \ ( a \ ( D i 



5.4. The Psuedo-Three-Body Resonance 

To complete the argument that the perturbations from the B ring and Mimas can act 
in concert to align the orbit pericenters in the Cassini Division, we will now demonstrate the 
existence of three-body-like resonances involving these objects. First, recall equation H4"l 
and realize that the overall goal is to find a term in d 2 w/dt 2 that is proportional to sin(ro — 
^B+jLt — <f ). To do this, we need to express d 2 w/dt 2 in terms of the disturbing function, 
which can be done by using the Lagrange planetary equations. 

Assuming the particles in the Cassini Division have no inclination, and keeping only 
terms to lowest order in eccentricity e. 

dw 1 dlZ 

17 = — 2~1~ + Wsec ' ( 17 > 
at na z e oe 

where w sec is again the precession induced by Saturn's oblateness, which will we assume 
here to be a constant at any given semi-major axis a. This means that: 

d 2 w d ( 1 dW 



dt 2 dt V na 2 e de 



(78) 
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In this case, we are only interested in the terms TZuimas and 1Zm ng derived above, and 
Tuning does not depend on the eccentricity e, so this equation simplifies to: 



d 2 W d ( 1 dlZMimas\ 



dt 2 dt \na 2 e de J 
Which can be expanded using Kepler's third law to give: 
d 2 w -1 



dt 2 na 2 e 



1 dTZuimas de 1 dTZuimas da d / dTZuimas \ 

e de dt 2a de dt dt \ de J 



(79) 



(80) 



Since de/dt and da/dt can be expressed in terms of derivatives of the disturbing function, 
the first two terms in this expression naturally give a product of two disturbing functions, 
enabling the term TZRi ng to mix with IZyiimas to obtain the desired frequency. However, 
the third term can also do this since: 

d ( dlZMimas\ d 2 TZMirnas de d 2 TZ Mimas da 

+ — snr: — H + 



dt \ de J dede dt dade dt 

Note that other derivatives are not of interest here because either they are zero or because 
they do not lead to terms with the appropriate frequencies. For example, d 2 lZMimas / {de) 2 = 
0, and any term involving dw /dt cannot give a useful term because TZmng is independent 
ofe. Thus: 



d 2 w -1 



dt 2 na 2 e 



1 dTZuimas de ( 1 dlZMimas d 2 TZMimas \ da d 2 1ZMimas de 



e de dt \2a de dade J dt dede dt 



(82) 



Now we can use the following Lagrange equations: 

da _ 2 dTZ 
dt na de 

de _ e dTZ 1 dTZ 
dt 2na 2 de na 2 e dw 



(83) 
(84) 



de 2 dTZ e dTZ 

dt na da 2na 2 de ' 
Note that we want terms involving a product of TZuimas and TZRi ng , so only derivatives of 
the latter will be of interest here. Furthermore, since TZjn ng only depends on A = nt + e 
and a, only the first term in each equation will contribute: 

Hence, substituting and simplifying: 



1 dTZuimas a d TZm imas 

\ dTZ Ring 2a 

imas 

dTZRing 



~dW ~ n 2 a*e [\2 de 2 dade J de dede da~\ ' 1 ' 

Notice that since TZuimas depends on various powers of a, d 2 TZMimas / dade is of order 
a~ l dTZMimas/de, so all the elements of the first term are of the same order. However, since 
TZm ng oc (a — Ro)~ x , dTZRing I 'da has terms that are of order (a — R )~ 1 TZRi n g, and since 
a — R is much less than a, this means the latter term is much larger, so we can approximate: 

d 2 W —2a d 2 TZuimas dTZRing 



dt 2 n 2 a A e dede da 



(87) 
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Substituting in Equations 1531 and 1691 into this expression (i.e., assuming = fj,2 = 0), 
we get: 



d 2 w 



2a 



Gm 



M 



/(a/ajw) sin(A + w — 2\mv 



2GaW a D 1 D 2 
3 R n (a — R n ) 



cos(A + wb 



Keeping only the desired long-period term, and recalling that 4>lm = 2\mi 
becomes: 



+ 2(j>L, this 



d 2 w 
~dW 



2n2 
3~e 



rriM a 

M S dM 



f(a/a M ) 



aWR Q a 2 aD x D 2 



sin(z*7 — wb + 2^>l) 



Note that for a 2:1 resonance aj 'om/(«/ / «m) ' 
1999), so this equation can be written as: 

d 2 tu 2 ~ 

-/ 8in(tu - ct7 B + 20zJ. 



M 5 R? (a-R f 

-0.75 (see Table 8.5 in Murray and Dermott 



where 



J o 



2n 2 
3~e 



dt 2 



rriM a 
M s a M 



\f(a/a M )\ 



aWR a 2 aD x t>i 
M s Rjia-Ro) 



(90) 
(91) 



is a positive quantity that indicates the strength of this resonance. 



Recall that compound trigonometric functions like sin(t%7 — wb + 20l) = sin(ro — wb + 
2(j>i J sm(Lt — Ol)) can be expressed as a series of terms of the form sin(ro — tub + jLt). We 
may therefore expand the right-hand side of Equation [90] into the series of terms: 



d 2 w 
~dW 



-fj sm(ru -w B + jLt). 



(92) 



where f 2 = Cjf 2 measures the strength of the individual resonances. As shown in Figure [TTT 
Cj will in general be a decreasing function of j, but for large- amplitude librations Cj would 
not be much less than unity for most j in the Cassini Division. 

If we consider the resonant argument ipj = w — wb +jLt, Equation [92] says that there 
is a term in the equation of motion: 



Since ipj 



dip i 



(ifj /2.), this means: 



-ff sin (fj. 



4/|sin 2 (^/2). 



(93) 



(94) 



where <p is the value of (fj where (fij=0. Note a librating solution to this equation of motion 
will only exist if (f>j = for some value of ifj. This will only happen if ip 2 < 4/?. Or, in 
other words, when: 

\w - w B - jL\ < 2fj (95) 

where w D is the apsidal precession rate when ip = 0. Now, the apsidal precession rate in 
the Cassini Division is dominated by Saturn's oblateness. Thus this resonance will only be 
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able to confine the pericenter locations over a region in the rings where the precession rate 
is within 2fj of wb — jL. Inserting the same numbers used in the previous section into 
Equation 1914 we find that in the simple case without mass anomalies (and the amplitude of 
the m = 2 pattern on the B ring edge is nearly constant): 



3 



f 2 ^7*10- 22 / S 2 cY^ ( - ) (-^-\ (^) (-M ( 10mkm ) 
Jj 1 J \aeJ \100g/cm 2 J \100km J \20km J J \ a - R Q J 

(96) 

so 2fj is of order 5* 10~ 11 /s or 0.00025°/day. Assuming the precession rates are dominated 
by Saturn's oblateness (so 5w/w ~ (7/2)5a/a), this range in precession rates corresponds 
to a range in semi-major axis on the order of a few kilometers. These resonances should 
therefore be able to align the pericenters of particle orbits over a region in the rings a 
few kilometers across, which is comparable to the amplitudes of the observed eccentric fea- 
tures. This result is encouraging, although we must point out that this analysis merely 
demonstrates that eccentric orbits within a few kilometers of the resonance can have their 
pericenters aligned by these resonances. This analysis does not determine the actual ec- 
centricity of individual particle orbits near these resonances, which must be the subject of 
future work. 



6. Discussion of Potential Further Work 

The above analysis shows that if the orientation of the B-ring edge librates at ~ 
0.06°/day, then the combined gravitational perturbations from the B ring and Mimas should 
align the pericenters of particle orbits at the observed locations of the eccentric gap edges 
in the Cassini Division. Furthermore, if the physical parameters of the B ring are close to 
the values assumed above, then particles with aligned orbital pericenters will extend over a 
range of semi-major axis comparable to the observed amplitudes of the eccentric features in 
the Cassini Division, which is needed to avoid streamline crossing. Both of these findings 
lend support to the notion that the gaps in the Cassini Division may indeed be generated 
by pseudo-three-body resonances involving perturbations from the Mimas and the B-ring 
edge. At the same time, it is also clear that more work needs to be done to develop a 
complete model of how the Cassini-Division gaps are formed in this sort of scenario. While 
such a complete model is beyond the scope of this paper, we will describe in this section 
what we see as productive avenues for future work towards this goal. 

First of all, while the above calculations provide an explanation for why there should 
be coherent eccentric features at the observed locations of the inner edges of the gaps in the 
Cassini Division, this analysis does not yet provide a physical explanation of why gaps should 
form at these locations. We expect that once the gaps have formed, the relevant resonances 
should maintain the eccentric shape of the inner edges by aligning the pericenters of the 
relevant particles' orbits. Furthermore, we speculate that the mass anomaly produced by 
the eccentric inner edge of each gap may help keep the outer edge from diffusing inwards via 
a modified shepherding mechanism in which the periodic radial and tangential gravitational 
forces from the eccentric inner edge cause the semi-major axes of the particles on the gap's 
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outer edge to migrate outwards. While iRappaportl ()1998l ) determined that the torque 



from a simple precessing eccentric ringlet is always negative and is therefore unable to 
confine the outer edge of a gap, the perturbations from eccentric edges (which could involve 
significant mass anomalies as well as well as radial excursions) have not yet been thoroughly 
investigated. Regardless of the specific shepherding mechanism, any successful model of the 
gap outer edges must be consistent with the observed near-circularity of the outer edges. 

While further analytical calculations could clarify which torques might be able to hold 
open these gaps, numerical studies may be needed to determine how the eccentric inner edges 
form in the first place. The above calculations suggest that the pericenters of eccentric orbits 
will become aligned, but they do not indicate whether there is any term in the equations 
of motion that would tend to generate eccentricities at these locations. While it is possible 
that such terms are produced by some combination of the various perturbations described 
above, it is also possible that the eccentricities are generated by an instability in the Cassini 
Division itself. If particle collisions produce small eccentricities, which are in turn aligned 
by the above resonance, this will produce small anomalies in the local gravitational field 
that could then drive larger eccentricities in the region, which will also become aligned, 
until an eccentric edge is formed. If we can better understand how these edges form, we 
may be able to convert the observed amplitudes of the radial excursions into constraints on 
the physical characteristics of the B ring such as its mass density. 

Whereas the amplitude of the edges' radial excursions should be related to the physical 
properties of the B-ring edge, the orientations of the m = 1 patterns should be related 
simply to the kinematics of the B-ring edge. In the idealized model given above, the aligned 
pericenter locations for the different edges are determined by the movements of the different 
components in the B-ring edge, and we expect this to be the case in the actual Cassini 
Division as well. However, the actual pericenter locations of the eccentric Cassini-Division 
features could be affected by finite mass anomalies in the outer B ring, as well as variations 
in the amplitude of the m = 2 pattern in the B ring edge. Given that these aspects of 
the B-ring edge are not yet well modeled, we cannot yet predict the absolute orientation 
of the Cassini-Division edges at any given time. In fact, if one accepts that the eccentric 
Cassini-Division edges are produced by the sorts of resonances described above, one could 
even use the observed amplitudes and orientations of the different eccentric edges to place 
constraints on the time variability in the B-ring's edge. 

More generally, the above model was developed assuming the m = 2 pattern in the 
B-ring edge librates with a frequency of around 0.06°/day, which matches the spacing 
in the observed pattern speeds of the eccentric features in the Cassini Division. More 
work needs to be done to confirm that this particular solution is indeed the correct one. 
Beyond considering more data (e.g. from Cassini images, Spitale and Porco in prep), more 
sophisticated theoretical modeling of the B-ring edge may be able to better constrain the 
likely solutions. Even if one of the other possible solutions for the motion of the m = 2 
pattern turns out to be correct, a variant of the above model may still be able to produce 
the observed pattern in the locations of the Cassini-Division features. For example, if the 
m = 2 pattern is actually circulating instead of librating, then its amplitude and drift rate 
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will periodically change as it drifts relative to Mimas, which should produce a series of 
resonant terms in the equations analogous to those derived above. 

Finally, a full understanding of the Cassini Division gaps will also need to encompass 
the noncircular ringlets in the Huygens, Herschel and Laplace gaps. The Huygens ringlet, 
being the most eccentric feature in the Cassini Division, may turn out to be the limiting case 
of an eccentric edge. The origin of the Lapla ce ringlet may be co nnected with the second- 
order Pandora resonance on its inner edge (jColwell et al.l l2009al ). Even so, the Herschel 
ringlet remains largely mysterious. 



7. Conclusions 

• The outer edges of five of the eight named gaps in the Cassini Division (that is, all 
the gaps that do not contain dense ringlets) are circular to within 1 km. 

• The inner edges of six of the eight gaps in the Cassini Division (that is, all of the gap 
inner edges that do not lie near first-order Lindblad resonances with known satellites) 
are eccentric, with m = 1 radial variations that drift around the planet at rates close 
to the expected apsidal precession rates. 

• The radial excursions in the inner edge of the Barnard Gap, which lies close to the 
5:4 ILR with Prometheus, do appear to have an m = 5 component tied to that moon. 

• The radial excursions in the outer edge of the B-ring, which lies close to the 2:1 ILR 
with Mimas, have an m = 2 component. However, the amplitude of the pattern and 
its orientation relative to Mimas change with time. While the available data do not 
yet determine a unique solution for the motion of the m = 2 pattern, one possibility 
is that it librates relative to Mimas with a frequency ~ 0.06°/day and an amplitude 
of ~ 50° in longitude (or a freqeuncy of ~ 0.12.°/day and an amplitude of ~ 30°). 

• In addition to the m = 2 pattern, the B-ring edge also appears to have an m = 
1 component that drifts around the planet at a rate close to the expected apsidal 
precession rate at the B-ring edge of 5.06° /day. 

• The pattern speeds of the eccentric features in the Cassini Division, including the 
Huygens ringlet and the m = 1 component of the B-ring edge, appear to form a 
regular series given by ft p = wb — jL, where j = 0, 1, 2, 3, 7, wb — 5.06°/day is 
the apsidal precession rate at the B-ring edge, and L ~ 0.06° /day is a likely value for 
the libration frequency of the m = 2 component in the B-ring edge. 

• By combining gravitational perturbations from both components of the B-ring edge 
with the perturbations from the 2:1 Mimas ILR, one can find terms in the equations 
of particle motion of the form zb oc sm(zv — tub + jLt). Such terms could act to align 
the pericenters of particle orbits at the locations of the eccentric inner edges of each of 
the gaps, and therefore could explain the locations of the gaps in the Cassini Division. 
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